Full counting statistics of chaotic cavities with many open channels 
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Explicit formulas are obtained for all moments and for all cumulants of the electric current through 
a quantum chaotic cavity attached to two ideal leads, thus providing the full counting statistics for 
this type of system. The approach is based on random matrix theory, and is valid in the limit when 
both leads have many open channels. For an arbitrary number of open channels we present the 
third cumulant and an example of non-linear statistics. 
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The physics of current fluctuations in mesoscopic con- 
ductors is an interesting and fundamental quantum me- 
chanical problem, since at low temperatures they are 
mainly due to the discreteness of the electron charge. 
The study of shot noise, for example, is an active area 
of theoretical and experimental research involving dif- 
ferent types of systems (e.g. quantum dots, disordered 
wires, quantum point contacts) and different regimes-'^ 
(e.g. Coulomb blockade, quantum Hall effect, localiza- 
tion). A relatively recent approach is the concept of full 
counting statistics^, the study of all cumulants of the 
charge fluctuation, which amounts to having complete 
information concerning charge counting in a transport 
process. This approach has recently attracted much at- 
tention and has been applied in a wide variety of situa- 
tions (see Ref. M and references therein) . Experimental 
measurements of the third moment of these fluctuations 
have already been reported^. 

In the case of chaotic cavities the random matrix the- 
ory (RMT) approach^ has been very successful in re- 
producing different experimental observations related to 
quantum transport, such as weak localization and uni- 
versal conductance fluctuations. If the cavity is con- 
nected to two ideal leads, supporting respectively A^i 
and N2 open channels, the conductance is given by the 
Landauer-Biittiker formula g — GoTr[T], where Go is 
the conductance quantum, T = tH and t is the trans- 
mission matrix. However, only very recently was an ex- 
act expression obtained^ within RMT for the shot noise 
P = PoTr[r(l-r)] (with Pq = 2eVGo where ^ is a smafl 
voltage bias) and an explicit general result is not avail- 
able for higher moments of the type Tr[T™]. For chaotic 
cavities with large A^i , N2 the third and fourth cumu- 
lants of charge transfer have been obtained'', as well as 
an expression for the cumulant generating function^. 

On the other hand, recent semiclassical calculations 
based on correlated classical trajectories that transmit 
through the cavity have been able to reproduce the RMT 
results, both for the conductance and for the shot noise. 
These calculations have a natural perturbative structure 
on the parameter N~^, where = iVi + A^2 is the to- 
tal number of channels. Initially the leading order ex- 
pressions were reproduced,— later the full series were 
obtained and exactly summed.— The next natural step 
would be to tackle Tr[T™], and it is thus of interest to 



have the corresponding RMT prediction of this quantity, 
at least to leading-order in N^^. This is the purpose of 
this work. 

We will be interested in the dimensionless moments 
defined as Y^^=i ^i™' where Ti are the eigenvalues of the 
matrix T and n = min{A^i, 7V2}- Within RMT the Ti are 
correlated random numbers between and 1 , whose dis- 
tribution depends only on the symmetries of the system 
(orthogonal, unitary or symplectic, labeled hy (3 — 1,2 
or 4 respectively). The average value of the moments are 
then given simply by 



M„ 



n{Tn. 



(1) 



The distribution of transmission eigenvalues can be 
characterized by a density, pp{T), such that (T™) = 

Pl3{T)T"^dT, or equivalently by a joint probability dis- 
tribution Vp such that 

(TD = dTi • • • dT„ TrVp{T). (2) 
Jo Jo 

The expression for V/3{T) is^ 



(3) 



where A(T) 
minant, a = 



Pi 



Y[i^j{Ti — Tj) is the Vandermonde deter- 
1 and Afg is a normal- 



2(|iV2 - iVil + 1) - i ana j\f3 
ization constant. In Ref. 6 the authors used simple re- 
currence relations from the theory of Selberg's integral— 
to obtain an exact result with arbitrary Ni , 7V2 for the 
second moment M2 and for the shot-noise (second cu- 
mulant). Here we follow a similar approach and com- 
pute the third cumulant, sometimes called the skewness. 
Moreover, we then proceed to obtain explicit formulas 
for all moments Mm and for all cumulants, valid to first 
order in the inverse number of channels, i.e. in the limit 
m,N2 » 1. 

We must note that in the semiclassical limit of short 
wavelengths some noiseless scattering states can be 
created, ^'^ leading to a breakdown of the universality 
implied by RMT predictions.— This phenomenon is 
governed by the ratio te/tjj of the quantum Ehrenfest 
time to the classical dwell time, and its influence has 



2 



been investigated on shot-noise, ^"'^^ the weak locahza- 
tion efFect^^ and conductance fluctuationsJ^^^^^^ Our re- 
sults are restricted to the universal regime te/td — > 0, 
when these system-specific corrections are neglected. 

Let us consider a certain fixed sequence of k positive 
integers, m — [mi, m^], and for any subsequence of 
length q < fc let us define the function 



(4) 



We take now TkP^{T)rfj{T), derive it with respect to 
Tfc and integrate over all variables to obtain 



F = {a + mk){P^{T)) + /? J] ( 
where the constant F is given by 



Tk 



(5) 



(6) 



We can see that F is actually independent of rrifc. 
Hence, we may equate the r.h.s. of ([5]) at different val- 
ues of this variable arriving at a recurrence relation. To 
solve this relation in general is presently beyond reach, 
but armed with some patience once can compute the first 
moments. This is essentially what was done in Ref.|^. We 
take it a bit further and find the third moment. Instead 
of writing the lengthy expression that arises for we 
present the corresponding cumulant (assuming for sim- 
plicity N2>Ni), 



Qs _ {N2-N, + l-j)iN2-N,-l + jj) 
Q2 (iV-l + |)(iV-3 + |) 

where Q2 is the average shot noise in units of Pq, 



(7) 



(P) 



iViA^2(iVi-l + |)(iV2-l + |) 



Po (vv_i + |)(7V-2+|)(7V~l + i)- 

The result for Q3 agrees in the limit ^ 1 with the one 
presented in Ref. 0- 

It is also possible to go beyond linear statistics, and 
compute higher correlations as for example 

n(n-i)(riT2(i-ri)(i-r2)) 
Q2 

_ (jVi~l)(jV2-l)(A^i-2+f)(iV2-2 + |) 

(iV-3 + |)(iV-4+|)(7V-2+|) ' ^> 

a quantity which would be important to compute the 
variance of the shot noise. 

To be able to arrive at a more general result, we now 
introduce the assumption that both leads contain a large 
number of open channels, A^i, N2 ^ 1, and thus n ^ k. 
In this case the main contribution to the summation in 



([5]) will come from j > k. We can thus approximate F 
by 



F^{a + mk){PL{T))+(3n(piiT) 



Tk 



(10) 



All the results obtained from now on should be under- 
stood as being valid to first order in N^^. Having said 
that, we drop the symbol and just write equalities. 
We can use the identity 



TV. 



1 /Tr-T„" 



(11) 



to simplify our expression for _F, 

F ^ ia + (3n){P^iT)}+p^ {P^-\T)Rt:.,{T)) , (12) 
where R^''^{T) denotes the symmetric polynomial 



(13) 



and we have neglected against a -\- (3n. 

Comparing (1121) for ruk and ruk — 1 we get the relations 

{PL-\T)n)=A2{P^-\T)}, (14) 

for TTifc ~ 1 and more generally 

{P^iT)) = {pL{T)T^') 

+ A, {p!^-'{T)[R^^I_,{T) - RtUiT)]) , 

(15) 

for mk > 2. In the previous equations Ai and A2 are the 
constants 

/3n Ni 2a + fin N2 

^1 = — -^TT = — - A2 = ^, , ^ , = — . (16) 



2(a + (in) N 



2{a + I3n) N 



Not surprisingly, the parameter (3 has dropped out of 
the calculation since leading-order results coincide for all 
universality classes. Iterating (jlSp k times we obtain 

(P^(r)) = A2{Pi^-\T)) A, {p^-\T)RtUiT) 

(17) 

Since we are interested in moments, we consider a par- 
ticular case of the previous equation which is 



(rn = A2-Ai(i?;'ii(r)) 

On the other hand, Eq. (fT7|) also gives 



(18) 



{Rl\T)) = A2Y,{Ti)-A, irr-'RfiT)) . (19) 

We must remark that the exponents of the terms inside 
the last brackets provide all ordered partitions of m into 
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feeding p3|) with our result ((22)) we can obtain the first 
few cumulants, and the coefficients -D^p are shown in Ta- 
ble II. We have found by direct inspection that these 
coefficients are such that 



QkiO = iV^(-l)'=+P^^^^^(fc - l,p - (24) 



where 



eral values of m and p. 



3 positive integers. These two equations can now be iter- 
ated together to yield the moments Mm = n{T"''), which 
will in fact be a polynomial of degree m, 



AL 



0=1 



NIN2 

7V2 ■ 



(20) 



Finding out the coefficient Cmp of the power is now 
an exercise in combinatorics. The first part of the prob- 
lem consists in answering the following question: In how 
many ways can one build sequences {ai, ...,a2p} with 
aj e {^1,^2} such that both Ai and A2 appear exactly 
p times and in all subsequences {oi, ...,a^}, q < 2p the 
number of A2^s is not larger than the number of Ais. 
The solution to this classic problem are the celebrated 
Catalan numbersr^ 



1 



'2p 
1 U 



(21) 



The power in Mm(^) will thus contain a factor 



(-1) 



■p— 1 



It will also be multiplied by another factor, 



which is equal to the number of ordered partitions of m 
into p positive integers. This is (™ri ) ■ 

We thus obtain our main result, an explicit expression 
for all the moments, valid to first order in the inverse 
number of channels: 



(22) 



The first three moments agree with known results.^ We 
present the coefficients Cmp with m up to 7 in Table I. 
The following equation. 



^ln{l + T,[e^-l]} = ^ 



fc=l 



k\ 



(23) 



relates the moments and the cumulantsi^ These will also 
be given by polynomials, Qfc(^) = ^^I]p=i ^fepC^- % 



(25) 



are the Stirling numbers of the second kind.^^ From the 
cumulants we can derive the generating function 



A 



fc=i 



dz 



1 + v/l + 4^(e-" - 1) 



(26) 



which is in fact equal to the one obtained in Ref. 0, thus 
implying the correctness of ([24]) .— 

In summary, we have explicitly obtained the random 
matrix theory prediction for all moments and all cumu- 
lants of the charge current in a chaotic cavity, in the 
limit of large channel numbers. Naturally, it would be 
desirable to obtain such explicit expressions for arbitrary 
channel numbers, but in this case we were able to com- 
pute only special cases such as (O and (0). The mo- 
ments are natural quantities to be studied in semiclassi- 
cal approaches to the problem, and indeed Eq. (|22|) has 
been reproduced using action-correlated trajectories in 
the open quantum star graph^^. The Hamiltonian case 
and corrections due to finite Ehrenfest time are discussed 
to some extent in Ref. H^- 
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TABLE II: The values of the cumulant coefficients Dtp for 
several values of k and p. 
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